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Abstract Vortices in a one-component dilute atomic ultracold Bose-Einstein con¬ 
densate (BEC) usually arise as a response to externally driven rotation. Apart 
from a few special situations, these vortices are singly quantized with unit circula¬ 
tion [5] . Recently, the NIST group has constructed a two-component BEC with a 
spin-orbit coupled Hamiltonian involving Pauli matrices [TU21I5] . and I here study 
the dynamics of a two-component vortex in such a spin-orbit coupled condensate. 
These spin-orbit coupled BECs use an applied magnetic field to split the hyperfine 
levels. Hence they rely on a focused laser beam to trap the atoms. In addition, two 
Raman laser beams create an effective (or synthetic) gauge potential. The result¬ 
ing spin-orbit Hamiltonian is discussed in some detail. The various laser beams are 
fixed in the laboratory, so that it is not feasible to nucleate a vortex by an applied 
rotation that would need to rotate all the laser beams and the magnetic field. In a 
one-component BEC, a vortex can also be created by a thermal quench, starting 
from the normal state and suddenly cooling deep into the condensed state [71. I 
propose that a similar method would work for a vortex in a spin-orbit coupled 
BEC. Such a vortex has two components, and each has its own circulation quan¬ 
tum number (typically 0, ±1). If both components have the same circulation, I find 
that the composite vortex should execute uniform precession, like that observed 
in a single-component BEC [73. In contrast, if one component has unit circulation 
and the other has zero circulation, then some fraction of the dynamical vortex 
trajectories should eventually leave the condensate, providing clear experimen¬ 
tal evidence for this unusual vortex structure. In the context of exciton-polariton 
condensates, such a vortex is known as a “half-quantum vortex” mm- 
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1 Introduction 

This study focuses on the possibility of creating and studying vortices in a two- 
component spin-orbit coupled Bose-Einstein condensate (BEC) of ultracold atoms. 
The principal motivation is the theoretical proposal by Spielman [1], and the sub¬ 
sequent two experimental papers by the NIST group: creation of vortices without 
rotation [2] and the study of spin-orbit structure in the Hamiltonian [3]- 

These experiments are modeled by a BEC in a thin essentially two-dimensional 
harmonic trap with tight confinement in the perpendicular direction. For a one- 
component condensate, many studies have shown the creation of vortices with a 
single unit of circulation, mostly by stirring the BEC to induce rotation [HE]. In 
a few cases, experiments have been able to study the real-time vortex dynamics 
for up to 1 s, typically a uniform precession m- 

The generalization of these ideas to a two-component spin-orbit coupled BEC 
involves the concept of synthetic gauge fields, and Sec. [2] contains an introduction 
to the NIST spin-orbit Hamiltonian. Here the principal emphasis is on a BEC in 
a trap provided by one or more focused laser beams. For completeness, this sec¬ 
tion also mentions the recent extensive studies of synthetic gauge fields in optical 
lattices (for both bosons and fermions), but they are not directly relevant to the 
present study of two-component vortices. 

Section [3] summarizes the experimental procedure [7] that used a sudden ther¬ 
mal quench from the normal cold atomic gas above the BEC transition temperature 
deep into the one-component BEC. Roughly 25% of the time, the quench created 
a singly quantized vortex (with random orientation), and the experiment was able 
to take 6-8 time-lapse images of the dynamical precession over nearly 1 s. 

As discussed in Sec. [5] the NIST procedure to create a spin-orbit coupled BEC 
uses a focused red-detuned laser beam to trap the cold atoms, a magnetic field 
along y to separate the F = 1 hyperfine levels, and a pair of opposing external 
Raman laser beams, typically aligned along ±x. These external laser beams and 
the magnetic field are fixed in the laboratory, which makes it difficult [8i to use 
the usual experimental approach of rotating the condensate to nucleate vortices 
in a one-component BEC (in a rotating frame, these external beams and fields 
become time dependent). In contrast, the thermal-quench method [7] should also 
work well to nucleate vortices in a two-component BEC because no rotation is 
required. 

Section [4] discusses my theoretical analysis that relies on the time-dependent 
variational Lagrangian formalism [9]. The two-dimensional vortex position 1*0 serves 
as a time-dependent variational parameter, and the resulting dynamical Lagrangian 
equations show that the vortex moves on a contour of constant energy. For a one- 
component vortex, this picture yields a uniform circular precession. In contrast, a 
two-component vortex offers a wider set of structures. Each component can have 
its own circulation quantum number mi and m 2 . If both have the same value (say 
+ 1), then the precession should remain uniform. Other possibilities exist, however, 
and a thermal quench should sometimes create a half-quantum vortex with (say) 
mi = 1 and m 2 = 0, which has been predicted and observed in excit on-polar it on 
condensates mm. My theoretical analysis suggests that a half-quantum vortex 
in a spin-orbit coupled condensate would have topologically distinct orbits, some 
of which would remain in the condensate, and others that would move to the edge 
of the condensate and then disappear. 
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2 Synthetic Gauge Fields 

Synthetic gauge fields and spin-orbit coupling have generated great excitement 
over the past decade especially with the initial realization by the NIST group m 
for a cold atomic gas of 87 Rb. This atom has a nuclear spin I = | and the single 
valence electron has a spin s = \ . Thus the vector sum F = I + s has two values 
F — 1 and 2. The lower manifold F = 1 has three substates mp = 1,0,—1. Of 
these, only the state |F, ttif) = j 1, —1) is confined in a typical magnetic trap. 

In practice, creating the effective spin-orbit coupling requires an applied mag¬ 
netic field to split the F = 1 manifold. Hence it is necessary to use an optical laser 
dipole trap (instead of a magnetic trap) to confine the cold atoms. For a given 
electric dipole moment d in an external electric field E, the energy is U = — d • E. 
A single neutral atom has an ac polarizability a(cj), which yields an induced dipole 
moment d(cj) = a(o;)E(a;). If the electric field is turned on adiabatically, the re¬ 
sulting energy is U = — ^a(uj)\E(uj)\ 2 (also known as the ac Stark effect). For 
low frequency (“red-detuned”) laser light, the polarizability is positive. Hence the 
atoms are drawn to regions of large \E\ 2 . In this case, a focused infrared laser 
beam will trap the atoms at the narrow waist, where \E\ 2 is largest. 

As an introduction to the idea of synthetic gauge fields, recall the transfor¬ 
mation to a frame rotating uniformly with an angular velocity ft. In this case, 
a simple analysis relates the Hamiltonian in the rotating frame H' to H in the 
laboratory (stationary) frame [5]: H' = H — ft • L, where L = r x p is the an¬ 
gular momentum. Rewrite the last term as — ft - L = — ft x r • p. Combine with 
free-particle kinetic energy to obtain 

_A_o xr n _ (P - MU x r) 2 M\n x r| 2 
2 M P 2 M 2 ^ 

where last term is a (negative) centrifugal potential that opposes any applied 
trapping potential Vt r . I can now interpret the term Mft x r as an effective gauge 
potential A e ff since it appears in the familiar combination (p — A eff ) 2 /(2M). More 
generally, whenever the Hamiltonian contains a term linear in p, the coefficient 
can be taken to define an effective (or synthetic) vector potential A e ff. 

In classical physics, a particle with charge q in a magnetic field B experiences 
a Lorentz force qw x B. For a quantum system, however, the focus is on the vector 
potential A, where B = V x A. In particular, when a charged particle moves from 
ri to r 2 along a path C, its wave function acquires a phase 

S=lJ r J A{v’)-dv'. (2) 

If it is possible to create such a phase, by whatever means, even a neutral par¬ 
ticle can experience a “synthetic” gauge field. Hence the new perspective is on 
“phase engineering” of the quantum state. Spielman at NIST has demonstrated 
synthetic gauge-field effects for trapped neutral atoms of 87 Rb m- More gener¬ 
ally, Refs. |121[T3lllT4lll5lll6j review this exciting and rapidly developing subject. 

Recently, various experimental groups have created synthetic gauge fields in 
optical lattices created by standing waves of interfering laser beams. In the tight- 
binding model for the lowest band, synthetic gauge fields appear as complex phases 
associated with the hopping parameters in the single-particle Hamiltonian. These 
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effects can arise by shaking an optical lattice with special time-dependent 
driving forces. This and related methods have created uniform synthetic flux in a 
two-dimensional optical lattice with 1/2 flux quantum per lattice plaquette [181119] . 
a realization of the topological Haldane model m in a distorted two-dimensional 
hexagonal optical lattice, and chiral edge states in Hall ribbons using “synthetic 
dimensions” [ 2022 ] • 

The creation of such synthetic gauge fields for atomic gases usually relies on 
strong laser fields that couple two or more atomic states. Typically, this coupling 
yields “dressed” eigenstates ^(r)), where the spatial dependence is crucial. When 
a particle moves adiabatically from n to r 2 , this spatial dependence yields a 
Berry’s phase 



(3) 


where A = ih(ip |Vt/j) is the synthetic gauge field. Here, we deal with neutral 
atoms, and it is convenient to take the effective charge as 1, so that A has the 
dimension of momentum. 

A typical normalized wave function has the form 




where both x and 77 depend on the spatial coordinate r. The resulting vector 
potential is 



(5) 


so that the state vector must have a spatially dependent phase 77 . Correspondingly, 
the induced synthetic magnetic field is B = V x A = — |^V(cosx) x V 77 . The 
essential conclusion is that we need ( 1 ) both 77 and x to have spatial dependence, 
and (2) their gradients must point in different directions. In practice, the NIST 
group [2] take 77 oc x and x ck y, leading to a synthetic “Landau” gauge, with 
A x oc y and a nearly uniform B Bz over a restricted region. 

The original paper [1 proposed a technique to create vortices in a non-rotating 
condensate, where the relevant angular momentum comes from the synthetic elec¬ 
tromagnetic field. The experiment [2] produced remarkable images of vortices 
shown as dark regions where the cores have reduced density. 

How does this example [2 work in detail? The discussion will lead to the 
important idea of spin-orbit coupling in an ultracold dilute atomic gas. A Bose- 
Einstein condensate is trapped in a red-detuned (typically infrared) focused laser. 
Apply two counter-propagating Raman laser beams along ±x, with ki,u;i and 
k 2 ,cj 2 , taking k\ « ~ ko (see Fig. [TJ . In a Raman transition, an atom absorbs 

a photon from one beam and emits a photon into the other beam, while making a 
transition between two different internal atomic states. The momentum transfer to 
the atom is « ±2hkox because of the recoil. Acoustic-optical modulators control 
the corresponding frequency transfer Auo = uji — 002 - 

In addition to the Raman laser beams along d=x, apply a Zeeman magnetic 
field along 7 /, splitting the three tuf states for the F — 1 manifold. It is possible 
to isolate the two states | 1 ) = | t) = 11, 0) and \2) = | = | 1 , — 1 ), which provides 

a convenient two-component basis. One can think of this pair as a pseudospin- \. 
In this two-component basis, the Raman beams act to couple the two pseudospin 
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Fig. 1 Condensate absorbs one photon with ki and emits another with k 2 from opposite 
Raman laser beams, thus acquiring a momentum ±2hkox. In addition, a Zeeman magnetic 
field along y splits the F = 1 manifold into substates, and our model focuses on only the two 
lower ones. Adapted from pQ with permission. 


states. Varying the applied magnetic field induces a detuning M/2 from the Raman 
resonance. 

In the present context, the most important effect of the Raman lasers is to 
induce an off-diagonal Rabi coupling between the two states | t) and | i). The 
relevant matrix element involves the electric dipole energy d • E, and the total 
electric field E has the spatial dependence e 2lk ° x , leading to the matrix element 


(t |d • E\ i) = ±hf2e 2ikoX 


( 6 ) 


Here the quantity i? is the Rabi frequency; it is fixed by the strength of the applied 
Raman laser beams. Note that I also neglect the effect of the trap potential and 
the Gaussian curvature of the trapping laser. 

In the two-component basis, the single-particle Hamiltonian becomes 



where I again omit the trap potential. The presence of spatially varying off-diagonal 
elements complicates the problem, but this spatial dependence can be removed 
with a unitary transformation m 



( 8 ) 


A simple analysis yields a new single-particle Hamiltonian Hso that now has a 
spin-orbit structure 


H so =U i H 0 U 


(9) 



( 10 ) 


| j^(7i 2 /M)(—iVT + koxcr z ) 2 + Ma z + Hf2a x ~^ 


(11) 
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where X denotes the 2x2 unit matrix and a J are the usual Pauli matrices. In effect, 
the Raman beams shift the minima of the two pseudospin dispersion relations to 
new and different local minima at = fhkox . These shifted minima represent the 
vector gauge fields (p — A) 2 , with A = —hkoa z x. In addition, the Rabi coupling 
induces the off-diagonal term ^hf2a x . 

To understand the new physics, continue to ignore the nonuniform trap poten¬ 
tial so that Hso has one-dimensional plane-wave solutions oc e lkx . Use /c^" 1 as the 
unit of length and the recoil energy Er = h 2 k^/2M as the unit of energy, leading 
to the dimensionless spin-orbit coupled single-particle Hamiltonian 


Hso 


f (k T l) 2 T 8/2 17/2 \ 

V n /2 (k-l) 2 -S/ 2 j' 


The associated eigenvalues follow immediately 

E±(k) = k 2 + 1 ± i v / (4fc + <5) 2 + r2 2 , 


( 12 ) 


(13) 


and I here focus on the lower band E-(k). 


2.1 limit of large Rabi frequency 


Two cases are of special interest, and the first is the behavior for large Rabi 
frequency (17 4). An expansion of E-(k ) in powers of 17 -1 yields 


E-(k) 


"2 


1 + 


17-4 

17 



overall shift effective mass gauge pote ntial shift 


(14) 


The first two terms are simply an overall downward energy shift, and the factor 
in front of the quadratic term is an effective mass. Note that the minimum in the 
dispersion relation is shifted from the usual position k = 0 to the new position 
6/(f2 — 4), identifying the x component of the synthetic vector potential as 


A 


X 


S 

17-4' 


(15) 


This is a central result of Spielman’s analysis [1 , namely the synthetic gauge 
field varies linearly with the detuning S. If 5 is constant, then there is no synthetic 
magnetic field because A x would be constant. To obtain a useful synthetic field, the 
experiment [2] uses a magnetic field gradient along y , so that 8{y) = S' y, with S' 
(a constant proportional to the field gradient) as a control parameter. In this way, 
the synthetic vector potential has the form of Landau gauge A x oc S' y, familiar 
from the quantum description of an electron in a uniform magnetic field. Its curl 
yields an effectively uniform synthetic magnetic field along z proportional to the 
control parameter S'. In this case, the neutral atoms experience an effective Lorentz 
force completely analogous to the real Coriolis force observed with a rotating 
condensate [Q. Note that we are here effectively already in the rotating frame, so 
that the vortices are at rest in the laboratory frame. 

Reference [2] used this approach to create vortices in a nonrotating condensate, 
as seen in Fig. [2] Note that these shapes become progressively more distorted with 
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Fig. 2 A synthetic vector potential A x = S' y in Landau gauge creates vortices for sufficiently 
large S'. Adapted from [2] with permission. 


increasing control parameter S ', which provides additional insight into the idea of 
a synthetic gauge field m- Here, we have A = —By x , where B is the synthetic 
uniform magnetic field. To record these images, the trap and (real) magnetic fields 
are suddenly turned off, which generates a synthetic electric field E = —dA/dt = 
Byx (ignoring any scalar potential). The resulting pulsed electric field produces 
an impulsive shear, as seen in Fig. [2] 


2.2 spin-orbit structure of Hso 


The same Hamiltonian in Eq. (11) exhibits rather different behavior for small 


values of the Rabi frequency i? < 4 Er (now in usual units), which emphasizes 
the spin-orbit character of the interaction. The cross term in the kinetic energy 
is linear in the momentum and has the form —p • A/M = p x hkoa z / M , which 
exhibits the matrix synthetic gauge field 


A x = —hko<J z 


(16) 


Note that this interaction is somewhat different from that familiar in atomic 
physics, which is proportional toL-cr = rxp-cr. It arises from similar structure in 
semiconductor physics, known variously as Rashba or Dresselhaus coupling. Since 
A here has only one component, there is no question of non-Abelian gauge fields, 
for that requires two or more noncommuting components of A. 

I now focus on the situation of zero detuning S = 0 and small dimensionless 


Rabi coupling, in which case the dimensionless Eq. (13) reduces to 


E±(k) =k 2 + 1 ± |^/l6fc 2 + fP. 


(17) 


If the Rabi coupling vanishes (i? = 0), this expression reduces to two shifted 
parabolas {k d= l) 2 that intersect at k = 0. For finite i2, however, an avoided 
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crossing splits the dispersion curves into an upper and a lower band. Reference |3] 
mapped out this behavior in their study of spin-orbit coupling in cold 87 Rb atoms, 
as shown in Fig. [3] Note the increasing splitting of the two bands with increasing 



- 2-1012 

Quasimomentum q/k L 


Fig. 3 Spin-orbit upper and lower bands for various values of 47. The gray intersecting parabo¬ 
las are for 47 = 0, and successive colored curves show the behavior with increasing 47. Adapted 
from ^ with permission. 


47. The two minima in the lower band become shallower and move closer together 
with increasing 47. For 47 < 4, the two local minima are at k 2 = 1 —47 2 /16, whereas 
for 47 > 4, there is only a single minimum at k = 0. Figure [4] verifies these results 
in great experimental detail [3j. 

This NIST scheme for creating synthetic gauge fields and spin-orbit coupling 
singles out x as a preferred drection. There are many proposals for more sym¬ 
metric spin-orbit coupling, and a typical case is pure Rashba coupling with the 
Hamiltonian 

Hr m km ( p x a y - p y a x ) /M, (18) 

where k is a coupling constant with the dimension of wave number. This form 
has two components of synthetic vector potential A x = Kkg v and A y = —hKa x . 
Note that [A x ,A y ] / 0 because the two Pauli matrices do not commute. In such 
a situation, these gauge fields are generally known as non-Abelian, and many 
unusual and intriguing properties can arise mmm- 

The lower band of the eigenvalue spectrum of Hr in Eq. ( |18| ) [including the 
free-particle term p 2 / (2M)T\ has a minimum on a circle of radius \p\ = fe, with a 
shape like the Mexican-hat potential. This special form has a linear crossing with a 
Dirac cone between the upper and lower bands, like an axisymmetric version of the 
gray curves in Fig. [3] In practice, the Rashba Hamiltonian can contain additional 
control terms, such as a detuning \HAcr z , which plays a role analogous to mass 
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Minima location in units of k L 

Fig. 4 Experimental trace of the position of the minima in the dispersion relation shown in 
Fig. [ 3 ] as a function of the Rabi frequency Q. Adapted from [3] with permission. 


in the Dirac theory and splits the upper and lower axisymmetric bands. Despite 
great experimental efforts, no such Rashba coupling has yet been achieved. 


3 Recent experiment on vortex dynamics in a single-component BEC 

For reasons that will become clear, it is valuable to review a recent experiment [7] 
that created vortices in a nonrotating condensate and took nondestructive time- 
lapse images of the subsequent vortex dynamics. A parabolic magnetic trap con¬ 
fined 87 Rb atoms in the particular hyperfine state 11, —1). The experiment started 
in the normal state and quenched rapidly into the superfluid state to low temper¬ 
ature T/T c < 0.4, with no measurable thermal cloud. 

Roughly 25% of the time, they found a vortex in the condensate with random 
(d=) orientation in a disk-shaped condensate. They applied a short microwave pulse 
that transferred « 5% of the atoms to the untrapped state 12, 0). These atoms fall 
under gravity and expand, allowing a direct image. Because these atoms were part 
of the original BEC, they provide a faithful small copy of the whole condensate, 
including the image of the vortex core. The experiment could repeat this process 6- 
8 times at intervals of ~ 90 ms, allowing a real-time study of the vortex dynamics. 
Figure [5] shows a dramatic set of images of a precessing vortex, with the upper 
row showing raw data and the second row the smoothed set of images. The lower 
part shows the fit to a uniform precession over approximately two full cycles. 

In addition to finding single vortices, they occasionally a few % of the 
quenches) found a vortex pair (often called a vortex dipole), which is a + vortex 
and a - vortex close together. For both the single vortex and the vortex pair, they 
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Fig. 5 A. Time-lapse images of a vortex taken every 90 ms, showing essentially two full cycles 
of precession around the disk-shaped trap. B. Angular position as a function of the time, 
showing highly uniform precession. Adapted from 7] with permission. 


found good agreement between the observed motion and that predicted with the 
Gross-Pitaevskii theory. 

Why is this particular experiment relevant for vortex dynamics in spin-orbit 
coupled condensates? A few years ago, Radic et al. [8] discussed vortices in such 
spin-orbit coupled systems and pointed out that it would be highly challenging to 
rotate not only the condensate but also the Raman laser beams and the magnetic 
field. At present, there have been no reports of such rotation experiments, and 
the thermal quench of Ref. [7 in principle provides a simple way to study vor¬ 
tices in such a nonrotating spin-orbit condensate. If this approach can indeed be 
implemented, it would also provide a detailed description of the associated vortex 
dynamics [9L 


4 Theory of vortex dynamics in a spin-orbit coupled Bose-Einstein 
condensate 

I rely on the time-dependent variational Lagrangian formalism that has proved 
valuable in studying the dynamics of vortices in trapped Bose-Einstein conden¬ 
sates HE]. Note that this variational analysis specifically includes both the confin¬ 
ing harmonic trap energy and the interaction energy, based on the Thomas-Fermi 
approximation. 


4.1 vortex dynamics in a one-component condensate 

Consider first a one-component cold atomic gas that is tightly confined in the z 
direction. In this case, it will form an effectively two-dimensional Bose-Einstein 
condensate with a condensate wave function ^(r). The Lagrangian here is given 
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by 


where 


L = T — E, 


T = — 
2 


j d \{. 


dt dt 


and E is the familiar Gross-Pitaevskii energy functional 

*2 


/* 


Em dr 


+ Vtrtf* + 


(19) 

( 20 ) 


( 21 ) 


Here, the three terms are the kinetic energy, the confinement energy of the trap, 
and the interaction energy with effective two-dimensional coupling constant # 2 D- 
Variation of the Lagrangian L with respect to E^ readily yields the exact time- 
dependent Gross-Pitaevskii (GP) equation. As usual with any variational principle, 
the current L also provides a valuable basis for a variational approximation. 

The strategy is to assume a trial wave function E(r, ro) with the two-dimensional 
position of the vortex ro as the time-dependent variational parameter. In partic¬ 
ular, I assume a normalized trial function 


E = 



( 22 ) 


where the first factor ensures the normalization f d 2 r\E\ 2 = N , and the second 
yields the Thomas-Fermi shape for the condensate density in a two-dimensional 
harmonic trap with condensate radius R. The position of the vortex ro = (xo,yo) 
appears in the phase 

Sir, ro) = arctan (——— ] . (23) 

y x xq J 


For the GP energy (21), ro affects only the kinetic energy, through the quantity 


VS = 


z x (r - rp) 
|r - r 0 | 2 


(24) 


which is essentially the induced flow velocity around the vortex core at ro. 

It is convenient to use plane polar coordinates ro = (ro, (j>o). A detailed calcula¬ 
tion yields the dimensionless energy E(uq), where uo = ro/R is the dimensionless 
scaled radial position of the vortex. For the one-component condensate, E does 
not depend on </>o. Similarly, the term T becomes 


T = -Nh4 >0 («0 - 14) • 


(25) 


The usual Lagrangian dynamics yields a pair of coupled equations 


1 


8E 


uo 


4> o = — 


uo (1 - Uq) d(j)o ' 
1 dE 


uo(l - Uq) duo 


(26) 


(27) 
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The Hamiltonian structure of these two equations ensures that dE/dt = 0 as the 
vortex executes its dynamical trajectory. Thus the vortex quite generally moves on 
a contour of fixed E. For a vortex in a one-component condensate, E is independent 
of 0o, so that iio = 0 and the motion is uniform circular precession with an angular 


velocity given by Eq. (27): 




(28) 


where £ is the vortex core radius. This result agrees well with experimental ob¬ 
servations [4]. Figure |6] shows typical contours of equal energy. A vortex precesses 
uniformly along such a circular curve with an angular velocity proportional to the 
radial gradient —dE/duo 



Fig. 6 Contours of constant dimensionless vortex energy E(uo) in a one-component conden¬ 
sate, where Uq = (x^ + y^)/R 2 is the vortex’s squared dimensionless radial position and R is 
the condensate radius. Adapted from [9] with permission. 


4.2 vortex dynamics in a spin-orbit coupled two-component condensate 


It is not difficult to generalize the Lagrangian to the more interesting case of a 
spin-orbit coupled condensate. Apart from the trap and interaction energies that 
remain unchanged, the new feature is the single-particle Hamiltonian of the form 
used in the NIST experiments [2,3] 


Eo 


(P - A) 2 

2 M 


hS z hfi x 

+ T < ’ + -T T 


(29) 


with A = — hkoxa z . Use of this Hamiltonian yields modified terms in the GP 
energy E+Eso = / d 2 r \P^'Ho\P. The experiment can control various parameters: 
the Raman laser wavenumber ko, the detuning 5 and the Rabi coupling strength 
12 . 

The trial function now has two components 


E = 




1/2 


c 


(30) 
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Here the first two factors are the same as in the one-component case (22), and £ 
is a two-component normalized spinor 


c 


e lSl cos(x/2) 
e lS2 e lr] sin(x/2) 


with two separate phases Si and S 2 [compare Eq. @] 

y-yo 


Sj = rrij arctan 


x — xq 


(31) 


(32) 


where rrij is an integer (typically rrij = 0, =bl). This structure assumes a vortex 
located at xo,yo, with quantized circulation mi, m 2 in the upper and lower com¬ 
ponents, respectively; in addition, the parameter a allows for an induced velocity 
along the preferred direction x. 

If mi s= m 2 with equal circulations, the resulting vortex dynamics is the same 
as for a single-component situation. If mi / m 2 (namely different circulations), 
however, the dynamics of the two-component vortex line is qualitatively different. 
I propose using a thermal quench like that of Ref. [7], anticipating that such an 
experiment would sometimes create a vortex with different circulations. 

The specific case mi = 1 and m 2 = 0 is analogous to a “half-quantum vortex” 
that has been predicted in thin films of superfluid 3 He-A and observed for both 
exciton-polariton BECs El and chiral p- wave superconductors m- A detailed 
analysis shows that the energy E now has terms proportional to cos0o and sin0o, 
meaning that the vortex dynamical trajectory now involves radial motion as well 
as azimuthal motion. Nevertheless, the vortex continues to move on a contour 
of constant energy E(uq,^>o). Unlike the previous situation, the vortex can have 
trajectories that leave the condensate along with those that remain inside. Figure [7] 
shows such contours for two different values of Rabi frequency O/Er = 0.1 (left) 
and i? /Er = 0.2 (right). 



Fig. 7 Contours of constant dimensionless vortex energy E(uo,0o) in a two-component con¬ 
densate with half-quantum vortex with circulations mi = 1 and m 2 = 0. Here, terms pro¬ 
portional to cos 0o and sin0o shift the center of the energy contours. The first figure shows 
contours for Q/Er = 0.1 and the second for Q/Er = 0.2. Adapted from [9] with permission. 
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5 Discussion and Conclusions 


This article arose from a presentation at a workshop on quantum gases, fluids, 
and solids, involving both the helium community and the cold-atom community. 
For that reason, I include a treatment of the recent cold-atom achievements in 
creating a two-component spin-orbit coupled Hamiltonian fni2ii3| . which will not 
be familiar to the broader low-temperature community. Owing to various lasers 
and magnetic fields that are fixed in the laboratory, it is not possible to transform 
to a single rotating frame with a time-independent Hamiltonian. Thus creating a 
two-component vortex must rely on other approaches, and I here propose a rapid 
thermal quench from the normal thermal cloud of cold atoms deep into the BEC. 
This method has successfully created singly quantized vortices in a one-component 
system |7], and the same technique should also work in a spin-orbit coupled BEC. 
In addition to conventional two-component vortices with the same circulation in 
both components, the two-component structure should also allow half-quantum 
vortices, in which one component has unit circulation and the other has zero 
circulation. If such half-quantum vortices exist, I find that their dynamics would 
be distinctive, in that a fraction of the vortex orbits would leave the condensate 
(see Fig. [7]). 

The variational approach in Sec. [4] has several inherent limitations. It assumes 
a Thomas-Fermi form for the density and spatially uniform spinor parameters x 
and r\. For a plane-wave solution, these parameters would depend on the wave 
vector k, and the present spinor yields the best constant variational choice. My 
wave function also takes the vortex singularity to have the same location for both 
components. In addition, the resulting vortex dynamics becomes singular near 
the outer edge of the condensate. To improve the description, a full numerical 
solution of the two-component Gross-Pitaevskii equation is probably preferable to 
a modified variational trial function. 


I have assumed small Rabi frequency Q/Er < 0.2 to ensure miscibility of the 
two components [3], but experiments for larger values would also be of interest. 
Since it seems necessary to use an optical trap, techniques are needed to release a 
small coherent fraction of the condensate atoms, but related methods have served 
well in similar contexts |24| . 


The NIST group created spin-orbit coupling in a trap, where the Raman beams 
provide a preferred direction. Nevertheless, many groups have proposed a more 
symmetric Rashba coupling, as seen in Eq. (18). It will be interesting to extend 
the current study to include the dynamics of vortices in the presence of such 
symmetric Rashba coupling. 


Acknowledgements Part of this article was written during a visit to the Institute for Ad¬ 
vanced Study, Tsinghua University, Beijing, and I am grateful to T.-L. Ho and H. Zhai for 
their hospitality. I thank W. Zheng for a valuable discussion concerning the interpretation of 
the synthetic electric field in Fig. [2] I. Spielman and D. Hall have provided copies of some of 
their figures, and I thank them for this assistance. I am grateful to G.-q. Liu and D. W. Snoke 
for discussions of the half-quantum vortices in exciton-polariton condensates. 




Vortex Dynamics in a Spin-Orbit Coupled Bose-Einstein Condensate 


15 


References 

1. I. B. Spielman, Raman processes and effective gauge potentials, Phys. Rev. A 79, 063613 
(2009). 

2. Y.-J. Lin, R. L. Compton, K. Jimenez-Garcia, J. V. Porto, and I. B. Spielman, Synthetic 
magnetic fields for untracold neutral atoms, Nature 462, 628-632 (2009). 

3. Y.-J. Lin, K. Jimenez-Garcia, and I. B. Spielman, Spin-orbit-coupled Bose-Einstein con¬ 
densates, Nature 471, 83-87 (2011). 

4. A. L. Fetter and A. A. Svidzinsky, Vortices in a trapped dilute Bose-Einstein condensate, 
J. Phys. Condens. Matter 13, R135-R194 (2001). 

5. A. L. Fetter, Rotating trapped Bose-Einstein condensates, Rev. Mod. Phys. 81, 647-691 
(2009). 

6 . B. P. Anderson, P. C. Haljan, C. E. Weiman, and E. A. Cornell, Vortex precession in Bose- 
Einstein condensates: Observations with filled and empty cores, Phys. Rev. Lett. 85, 2857 
( 2000 ). 

7. D. V. Freilich, D. M. Bianchi, A. M. Kaufman, T. K. Langin, and D. S. Hall, Real-time 
dynamics of single vortex lines and vortex dipoles in a Bose-Einstein condensate, Science 
329, 1182-1185 (2010). 

8 . J. Radic, T. A. Sedrakyan, I. B. Spielman, and V. Galitski, Vortices in spin-orbit coupled 
Bose-Einstein condensates, Phys. Rev. A 84, 063604 (2011). 

9. A. L. Fetter, Vortex dynamics in spin-orbit coupled Bose-Einstein condensates, Phys. Rev. 
A 89, 023629 (2014). 

10. Y. G. Rubo, Half vortices in exciton polariton condensates, Phys. Rev. Lett. 99, 106401 
(2007). 

11 . K. G. Lagoudakis, T. Ostatnicky, A. V. Kavokin, Y. G. Rubo, R. Andre, and B. Deveaud- 
Pledran, Observation of half-quantum vortices in an exciton-polariton condensate, Science 
326, 974-976 (2009). 

12 . J. Dalibard, F. Gerbier, G. Juzeliunas, and P. Ohberg, Colloquium: Artificial gauge po¬ 
tentials for neutral atoms, Rev. Mod. Phys. 83, 1523-1543 (2011). 

13. H. Zhai, Spin-orbit coupled quantum gases, Int. J. Mod. Phys. B 26, 1230001 (2012). 

14. V. Galitski and I. B. Spielman, Spin-orbit coupling in quantum gases, Nature 494, 49-54 
(2013). 

15. N. Goldman, G. Juzeliunas, P. Ohberg, and I. B. Spielman, Light-induced gauge fields for 
ultracold atoms, arXiv: 1308.6533v2 (2013). 

16. H. Zhai, Degenerate quantum gases with spin-orbit coupling, arXiv: 1403.8021vl (2014). 

17. J. Struck, C. Olschlager, M. Weinberg, P. Hauke, J. Simonet, A. Eckardt, M. Lewen- 
stein, K. Sengstock, and P. Windpassinger, Tunable gauge potential for neutral and spinless 
particles in driven lattices, Phys. Rev. Lett. 108, 225304 (2012). 

18. M. Aidelsburger, M. Atala, M. Lohse, J. T. Barreiro, B. Paredes, and I. Bloch, Realization 
of the Hofstdter Hamiltonian with ultracold atoms in optical lattices, Phys. Rev. Lett. Ill, 
185301 (2013). 

19. H. Miyake, G. A. Siviloglou, C. J. Kennedy, W. C. Burton, and W. Ketterle, Realizing the 
Harper Hamiltonian with laser-assisted tunneling in optical lattices, Phys. Rev. Lett. Ill, 
185302 (2013). 

20. G. Jotzu, M. Messer, R. Desbuquois, M. Lebrat, T. Uehlinger, D. Greif, and T. Esslinger, 
Experimental realization of the topological Haldane model, Nature 515, 237-240 (2014). 

21 . M. Mancini, G. Pagano, G. Cappellini, L. Livi, M. Rider, J. Catani, C. Sias, P. Zoller, 
M. Inguscio, M. Dalmonte, and L. Fallani, Observation of chiral edge states with neutral 
fermions in synthetic Hall ribbons, arXiv: 1502.02495vl (2015). 

22. B. K. Stuhl, H.-I. Lu, L. M. Aycock, D. Genkina, and I. B. Spielman, Visualizing edge 
states with an atomic Bose gas in the quantum Hall regime, arXiv: 1502.02496vl (2015). 

23. J. Jang, D. G Ferguson, V. Vakaryuk, R. Budakian, S. B. Chung, P. M. Goldbart, and Y. 
Maeno, Observation of half-height magnetization steps in Sr 2 Ru 04 , Science 331, 186-188 
( 2011 ). 

24. A. Ramanathan, S. R. Muniz, K. C. Wright, R. P. Anderson, W. D. Phillips, K. Helmerson, 
and G. C. Campbell, Partial-transfer absorption imaging: A versatile technique for optimal 
imaging of ultracold gases, Rev. Sci. Instrum. 83, 083119 (2012). 



